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Abstract 

We simulate 3-flavour lattice QCD at small isospin chemical potential fij and finite temperature 
T. At uj = there is a critical mass m c where the finite-temperature transition changes from 
first order to a crossover. We measure the /i/ dependence of the transition (3 (/3 C ) for m close to 
m c . (3 C and hence T c decrease slowly with increasing [ij. f3 c at finite \ij is in good agreement (5 C at 
finite fi (quark- number chemical potential). We use fourth-order Binder cumulants to determine 
the nature of this transition and to search for a critical endpoint. We measure the dt dependence 
of these cumulants and extrapolate to dt = 0. (dt is the 'time' increment used in the hybrid 
molecular-dynamics simulations.) Preliminary measurements of these Binder cumulants show lit- 
tle uj dependence. (Simulations at imaginary u indicate that the u dependence of the Binder 
cumulants is also weak.) This contrasts to the uj dependence we observed at fixed dt. 



* Talk presented by J. B. Kogut at the Workshop on QCD in Extreme Environments, Argonne National 
Laboratory, 29th June to 3rd July, 2004 



I. INTRODUCTION 



Although this is nominally the writeup of the talk delivered by JBK at the Workshop on 
QCD in Extreme Environments, we have taken the liberty of updating our results to reflect 
the 'data' from simulations performed since then. We have done this, since studies of the 
dt dependence of observables indicate that the Binder cumulants for the chiral condensate 
have a strong enough dependence on dt to change our conclusions. For our earlier results 
one should see our contribution to Lattice 2004, Fermilab jjj. 

QCD at finite baryon/quark- number density and temperature is relevant to the physics 
of relativistic heavy-ion collisions. The various accelerators which probe this region of the 
phase diagram of nuclear matter include RHIC, CERN heavy-ion, AGS and SIS. Between 
them they should map out a considerable region of this phase diagram. 

QCD and in particular Lattice QCD, at a finite chemical potential /z for quark-number, 
have a complex fermion determinant with a real part of indefinite sign. For this reason stan- 
dard simulation methods based on importance sampling cannot be applied directly. Various 
methods have been developed to circumvent this "problem for small u in the neighborhood 
of t fie finite t e m pe rature tr a nsltl o n ESflflflfiflflQ.Tfie transltlon * and We 
temperatures for finite /i and for finite \ij appear to be the same provided /// = 2/i. This 
equality was noted by the Bielefeld- Swansea group \^ and observed by us for the 2-flavor 
case Thus the phase of the determinant appears to be unimportant for determining 

the position of the transition. We have argued that this was because the fluctuations in the 
phase of the fermion determinant are moderate for small /i and high temperatures, on the 
lattice sizes we use. We speculate that the nature of the transition might also be insensitive 
to the phase of the determinant. If so, simulations at finite fii offer an alternative way to 
study the finite temperature transition at finite /i. Lattice QCD at finite fii has a positive 
fermion determinant, so standard simulation methods work. 

We are now extending our simulations to the 3-flavor case. Here we ignore the compli- 
cation that at finite /i/, this theory really describes QCD with 3/2 u quark flavors and 3/2 
d quark flavors. This theory can be interpreted as 3-flavor QCD at finite /i = Hi/2 where 
we use only the magnitude of the fermion determinant and ignore the phase. 3 flavors are 
of interest because they are closer to the physical 2 + 1 flavors. In addition, for quark mass 
m < m c , the phase transition is first order at \i — \i% — becoming a crossover for m > m c 
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I2I Q| . It is argued that the critical point at m = m c moves continuously to higher 
quark mass as \x is increased, where it becomes the sought-after critical endpoint p, si- 
Choosing a quark mass just above m c , one should be able to keep the critical endpoint as 
close to fi = as desired. In particular, one should be able to keep the critical \x in the 
range where the [i and fij transitions are related. Thus we might expect a critical endpoint 
at small /i/, which we could determine by direct simulations. 

Our initial simulations appeared to show evidence for this critical endpoint How- 
ever, the Binder cumulants, which we use to determine the nature of the transition, depend 
strongly on dt. We are currently performing simulations to study this dependence to allow 
extrapolation to dt = 0. Indications are that the fij dependence of the dt = cumulants 
is rather weak. This observation is consistent with what was predicted from analytic con- 
tinuation from imaginary /i by de Forcrand and Philipsen sj]. However, our preliminary 
measurments favor a slow increase with increasing fij while their measurements favor a slow 
decrease with increasing /i. Considerable more work will be needed to identify the critical 
endpoint if it exists. 

In section 2 we introduce our formulation of lattice QCD at finite \xi and discuss its 
relationship to QCD at finite \x. Section 3 describes our current simulations of 3-flavor 
lattice QCD at finite fii and gives preliminary results. In section 4 we give our conclusions 
and discuss what remains to be done. 



II. QCD AT FINITE m 

The staggered quark action for lattice QCD at finite isospin chemical potential /ij is 
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S f = X \x[lp(\ T ^i) + m }x + iAe^x] , (1) 

sites ^ 1 ' 

where p^^r^fij) is the normal staggered pfor a staggered fermion field x with 2 explicit flavor 
components, and r are the isospin Pauli spin matrices acting on this flavor space, with the 
links in the +t direction multiplied by exp(|r3yU/) and those in the —t direction multiplied 
by exp(— |t3//j). Explicit I3 symmetry breaking is provided by the term proportional to 
A. A 7^ is only needed for fij > m n where, in the low temperature phase, ^3 is broken 
spontaneously by a charged pion condensate. 
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The fermion determinant is 



det \[p^m) + m^lp^m) +m} + A 2 } 



(2) 



where in this case p is only a 1 x 1 matrix in isospin space, which means that we only need 
use a single flavour-component fermion field in our simulations. This determinant is real 
and positive allowing us to use standard hybrid molecular-dynamics simulations with noisy 
fermions, to allow us to tune the number of flavors from 8 down to Nf. Here we will choose 
Nf — 3, even though this represents 3/2 isodoublets and probably does not have a sensible 
continuum limit. 

When A = 0, this determinant is the magnitude of the fermion determinant at finite 
fj, = so simulating at finite /x/ is equivalent to simulating at finite /x using only the 
magnitude of the determinant and neglecting the phase. Hence, the theory we are considering 
corresponds to the 3-flavor theory at quark number chemical potential /x, in which the phase 
of the fermion determinant has been set to zero. We will argue that for small enough /x, the 
phase 6 of the determinant is well enough behaved, that it does not affect the position of 
the finite temperature transition. 

For simplicity let us consider a gluonic observable X. Its expectation value at finite /x is 
given by the ratio 



is reasonably large and varies smoothly as (3 crosses the transition at fixed /xj, i.e. it varies 
considerably more slowly than (X) for the volume considered, it will not have much effect 
on the transition. We can therefore consider e ld X as our observable and, to the extent to 
which the position of the transition is independent of the observable, the finite /x transition 
and the finite /xj transition will be coincident. It is not unreasonable to assume that the 
nature of the two transitions will also be the same. Because this argument only works for 
finite volumes, it should be considered as a plausibility argument rather than a proof. 




(3) 



where 6 is the phase of the fermion determinant. If the denominator 



(cos0)U =2/1 



(4) 
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III. SIMULATIONS AT SMALL m AND FINITE TEMPERATURE 



We are simulating the 3-flavor theory at finite temperature and small /i/, i.e. /i/ < m W} 
since the transition to the phase where isospin (7 3 ) is broken by a charged pion condensate 
occurs at jii = m n for zero temperature. Since no such condensed phase exists for QCD 
at finite fi, the physics at finite /i and that at finite /!/ can at best be related for /!/ < m n 
(/i < m,/2). 

Our simulations are being performed on N t = 4 lattices where the zero-chemical-potential 
finite-temperature phase transition is first order for m < m c (0) and a crossover for m > 
m c (0), with m c (0) ~ 0.033. At m = m c (0), the transition is a critical point in the universality 
class of the 3-dimensional Ising model. As fi is increased, m c is expected to increase so that 
the m c (/i) becomes the critical endpoint. Fixing m to lie just above m c (0), the transition 
should be a crossover for // < \i c with /i c (m) the solution of m = m c {jj, c ). For jj, > /i c , 
the transition should become first order. The critical endpoint should again lie in the Ising 
universality class. From the observation that the /i and fii transitions seem to be related, 
one might expect a corresponding endpoint for finite /i/. 

We are performing simulations at m = 0.03, m = 0.035 and m = 0.04. Most of our 
effort has been concentrated on m = 0.035, just above m c (0), where we might hope to find 
a critical endpoint for small /xj. Our simulations are being performed on 8 3 x 4, 12 3 x 4 and 
16 3 x 4 lattices. 

Figure Q shows the Wilson line from our 8 3 x 4 simulations at m = 0.035 for several 
values of /i/ < m n . The values of dt used for each of the [ij values used in these simulations 
are dt = 0.05 for ^ = 0.0, 0.2, dt = 0.04 at /i = 0.25, dt = 0.03125 for ^ = 0.03 and 
dt = 0.02 at [ii = 0.0375. The length of the simulations for each of those (3 values closest to 
the transition for each fij is 160,000 time-units(trajectories). 

We study the nature of these transitions using 4-th order Binder cumulants Q|. If X is 
an observable, its Binder cumulant is defined by 

B (X) - « X ~< X » 4 > (5) 

The Binder cumulant of interest is that for the chiral condensate. Since we only measure 
stochastic estimators of the chiral condensate, we need at least 4 independent stochastic 
estimates of (ipip) for each configuration to obtain an unbiased estimate of B4 - we actually 
use 5 noise vectors for each configuration. 



SU(3) N f =3 m=0.035 A=0 8 3 x4 lattice 



i — i 
O 

m 




P 



FIG. 1: Wilson Lines as functions of j3 for several values of on an 8 3 x 4 lattice with m = 0.035. 



On a lattice of infinite spatial extent, B^tpip)) would be 3 for a crossover and 1 for a 
first order transition, while for a critical point in the Ising universality class it would have 
the value 1.604(1). For finite lattices B^{{ipip)) will lie above the Ising value in the crossover 
region and below this value in the first-order regime, passing smoothly through a value close 
to the Ising value at the critical endpoint. As the size of the spatial lattice is increased, the 
curve will pass more steeply through the Ising value. The critical endpoint can be obtained 
by observing where the curves for different lattice sizes cross, or where the curve for a given 
large lattice passes through the Ising value. Ideally one should use several spatial volumes, 
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observe where the curves cross, and check that the value of B4 at this point is (close to) the 
Ising value. 

Preliminary indications based on simulations with the dt values given above showed some 
evidence for such a critical endpoint. However, it was observed that, while the value of B4 fell 
with increasing \ii at /// << m n , it began to increase again long before m n . This suggested 
that we should study the dt dependence of B4. We are now performing a detailed finite dt 
analysis. What we find is that the dt dependence of the fluctuation quantities - B^^ip)) 
and the susceptibility x({i J 4')) ~ is considerable, something we had also observed in the 2 
flavor case at Hi > m n . The preliminary results of this analysis which we present below 
indicate that considerable work will be needed to determine whether there is evidence for a 
critical endpoint, and if so, determine its position. 

Figure 121 shows the dt dependence of B± for the chiral condensate for m = 0.035 and 
Hi = from our 12 3 x 4 runs at dt values of 0.0625, 0.05, 0.04 and 0.03125. Each point is 
the value of B4 obtained from performing a Ferrenberg-Swendsen extrapolation 13] to its 
minimum from f3 values close to the transition. We have noted that the values of (3 C obtained 
from the minima of B4 and from the maxima of the chiral susceptibility are consistent. We 
have averaged the results obtained from 2 or 3 different (3 values close to the transition, for 
each dt. The length of the runs at each (3 was 160,000 At = 1 trajectories. We also noted 
that the values of /3 C obtained from different j3s and the same dt were consistent. We note 
considerable variation of B4 over the range of dts used in these runs. A fit linear in dt 2 
gives B A (dt = 0) = 1.84(3) with a x 2 /dof of 0.6. This shows that, for dt = 0.0625 the finite 
dt error « 16.5%. Even at dt = 0.05 which we used for our earlier simulations the error 
~ 9%. This also shows that finite dt errors can take one from B4 > 1.604, which is normally 
interpreted as a crossover to B4 < 1.604 which is commonly interpreted as indicative of a 
first order transition. 

We are performing similar calculations of the finite dt dependence of B4 for the chiral 
condensate for fii > 0. Preliminary results are shown in figure El Again we see a strong 
dependence of B4 on dt. Because hi 7^ makes the Dirac operator more singular, smaller 
dt values are required as Hi is increased. An extrapolation linear in dt 2 to dt = gives 
B4 = 1.89(3) for Hi = 0.2 and _B 4 = 2.01(8) for hi — 0.3. The rise in B4 with increasing fii, 
although suggestive, is not statistically significant. All we can claim at present is that there 
is little if any dependence of B^iptp)) on Hi over this range. For completeness we plot these 
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SU(3) 12 3 x4 lattice N f =3 m=0.035 /ll=0 X=0 
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FIG. 2: Estimates of B^-ipip)) as a function of dt 2 at m = 0.035 and m = on a 12 3 x 4 lattice. 

extrapolated values of in figure EJ 

The values of /3 C , the transition j3 are calculated from the minima of B±. We have noted 
that these estimate of (3 C are in excellent agreement with those obtained from the maxima 
of the susceptibilities for the chiral condensate. Extrapolating these to dt = 0, linearly in 
dt 2 , we obtain the results shown in figure The statistical errors are small, but there are 
large systematics associated with this preliminary 'data'. The straight line fit shown in the 
figure is: 

P c = 5.15195 -0.19428/4. (6) 

We do not quote any errors because of the preliminary nature of this result. If we interpret 
/ij = 2/x as we have suggested, this is in good agreement with the measurements of de 
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SU(3) 12 3 x4 lattice N f =3 m=0.035 ^,=0.2 X=0 



il 

pq 



m 




0.0000 0.0005 0.0010 0.0015 0.0020 0.0025 0.0030 

dt 2 

SU(3) 12 3 x4 lattice N f =3 m=0.035 ^,=0.3 X=0 
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FIG. 3: Estimates of B^iptp)) as a function of dt 2 at m = 0.035 on a 12 3 x 4 lattice a) for \ij = 0.2 
and b) for m = 0.3. 



Forcrand and Philipsen 

It is helpful to consider why the transition appears more like a first-order transition as dt 
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FIG. 4: Estimates of B^ipip)) extrapolated to dt = 0, as functions of fij on a 12 3 x 4 lattice. 

is increased from zero. A brief discussion of this was given in our 2-flavor paper 11]. One of 
the effects of simulating at a finite dt using the hybrid molecular-dynamics method, which is 
an inexact algorithm, is that the effective (3 value (3 e ffecUve measured from the kinetic term 
in the molecular-dynamics Hamiltonian using the equipartition theorem is (on the average) 
less than the input (3, so that we are effectively running at a lower (3 value. The quantity 
Aft = (3 — I3 e ff ec uve increases with increasing dt. This would not be a problem if the shift 
were roughly constant as (3 was varied across the transition. However, this is not the case 
for a constant dt. A(3 is significantly larger below the transition than above. We illustrate 
this from our simulations at /ii — 0, m — 0.035, A = and dt = 0.05 on an 8 3 x 4 lattice, 
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12 3 x4 N f =3 m=0.035 A=0, /8s at minima of B 4 
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FIG. 5: p c extrapolated to dt = as a function of /// on a 12 3 x 4 lattice. 

where we have 'data' at a number of f3 values over a range which includes the transition. 
This is shown in figure |HJ 

Now let us consider what the effect of such f3 dependence for A/3 has on the behavior 
of the simulation close to the transition. Here, a small change in j3 causes a relatively 
large change in (^effective- This causes a relatively large change in all observables. This, 
in turn, makes the transition appear more abrupt and hence nearer to being first-order. 
Looking at this in another way, ^effective needs to undergo a relatively large fluctuation to 
take it from one side of the transition to the other. This makes it more difficult for the 
system to fluctuate between its low-/3 domain and its high-/3 domain. Such a tendency to 
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8 d x4 lattice /a=0 m=0.035 X=0 dt=0.05 
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FIG. 6: A/3 = — f3 e ffective as a function of /3 at /// = 0, m = 0.035, A = and dt = 0.05 on an 
8 3 x 4 lattice. 



remain in one domain or the other, only 'tunneling' between domains rather infrequently, is 
a property one associates with first-order behavior. Hence the transition will appear to be 
closer to a first-order transition than it really is. These considerations are in accord with 
our observations. 
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IV. DISCUSSION AND CONCLUSIONS 



We are simulating 3-flavor lattice QCD with a chemical potential /// for isospin (J 3 ), at 
small Hi and temperatures close to the finite temperature transition from hadronic/nuclear 
matter to a quark-gluon plasma. These simulations are being performed on 8 3 x 4, 12 3 x 4, 
and 16 3 x 4 lattices. In this regime, the behavior of QCD at finite fii and that at finite 
quark-number chemical potential fi appear similar for fij = 2/i. QCD at finite \x and QCD 
at finite fii are related if one can neglect the phase angle of the determinant, or include it in 
the observables. The decrease in (3 C with increasing \ij that we have observed is consistent 
with that seen by de Forcrand and Philipsen if we make the substitution /i = which 
adds further support to the proposition that finite ji and finite fij are related. 

At zero chemical potentials there exists a critical mass m c . For m < m c the finite 
temperature transition is first order. For m > m c this transition is a mere crossover, with 
no real phase transition. At m = m c the transition is a critical point in the universality class 
of the 3-dimensional Ising model. For staggered lattices with N t = 4, m c ps 0.033. It has 
been suggested that at finite /i, the critical mass m c (/i) moves continuously (in jj) to higher 
values, where it becomes the critical endpoint. Hence the best place to look for a critical 
endpoint in ji or fij should be for a quark mass just above m c where it would be expected to 
be found at a small chemical potential. For this reason most of our simulations have been 
performed at m = 0.035. 

We have used 4-th order Binder cumulants for the chiral condensate to determine the 
nature of the finite temperature transition. A critical endpoint would be where the Binder 
cumulants for different lattice sizes cross. At this point its value should be close to the 
Ising value 1.604(1). Although our preliminary results appeared to show evidence for such 
an endpoint, further simulations studying the finite dt dependence of this cumulant showed 
the finite dt errors to be large at the dts we were using. Our current estimates show the 
extrapolated cumulants to lie well above the Ising value for the \ii range being considered. 
There is no indication that these cumulants are falling with increasing /i/. In fact our 
measurements show a slight rise with increasing However, this rise is not statistically 
significant, especially considering that these preliminary measurements could have sizable 
systematic errors. All that we can claim is that the Hi dependence of the Binder cumulants 
is very weak over this range. This is in qualitative agreement with what was observed at 
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finite \i by de Forcrand and Philipsen, who observed a slight decrease in B± with increasing 
/i 2 , but could not exclude a which is independent of fi 2 . 

If we do not find evidence for a B^ which decreases with increasing /jj at this quark mass, 
we should check the situation even closer to m c . The lack of any critical end point, even for 
m much closer to (but still above) m c would suggest that similar behavior would be found 
at small /z, since our estimates of the fluctuations of the phase angle as determined from 
(cos 6) {6 is the phase of the fermion determinant at \x = indicate that it is reasonably 
well behaved for m = 0.035 over most of the region fij < for an 8 3 x 4 lattice. This, in 
turn, would suggest that the critical endpoint at finite /i (if it exists) is fluctuation induced. 
Thus we would expect to find it where fluctuations in 9 become large, i.e. for /i > 
If so, it would not move continuously away from fi = as m becomes larger than m c , but 
might be expected to jump from \i — to \i ~ as m is increased beyond m c . 

The strong dependence of B^ on dt suggests that we should investigate the use of recently 
formulated exact algorithms (see and references contained therein), for our simulations. 
Since such methods rely on stochastic approximations to the fermion determinant, they 
could also help us in applying reweighting by the phase of the fermion determinant as a way 
of studying finite [i in a quantitative fashion, to obtain the equation-of-state. Here we note 
recent work indicating that finite fii simulations provide a better platform for reweighting 
than reweighting from simulations at \i — 
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